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ABSTRACT 

The merger of binary neutron stars is likely to lead to differentially rotating remnants. In this paper, we 
survey several cold nuclear equations of state (EOSs) and numerically construct models of differentially 
rotating neutron stars in general relativity. For each EOS we tabulate maximum allowed masses as 
a function of the degree of differential rotation. We also determine effective polytropic indices, and 
compare the maximum allowed masses with those for the corresponding polytropes. We consistently 
find larger mass increases for the polytropes, but even for the nuclear EOSs we typically find maximum 
masses 50% higher than the corresponding values for nonrotating (TOV) stars. We evaluate our findings 
for the six observed binary neutron star (pulsar) systems, including the recently discovered binary pulsar 
J0737-3039. For each EOS we determine whether their merger will automatically lead to prompt collapse 
to a black hole, or whether the remnant can be supported against collapse by uniform rotation (possibly 
as a supramassive star) or differential rotation (possibly as a hypermassive star). For hypermassive stars, 
delayed collapse to a black hole is likely. For the most recent EOSs we survey the merger remnants can all 
be supported by rotation against prompt collapse, but their actual fate will depend on the nonequilibrium 
dynamics of the coalescence event. Gravitational wave observations of coalescing binary neutron stars 
may be able to distinguish these outcomes - no, delayed or prompt collapse - and thereby constrain 
possible EOSs. 

Subject headings: Gravitation — relativity — stars: rotation 



1. INTRODUCTION 

The remnant formed in the coalescence of binary neu- 
tron stars is likely to be differentially rotating (see, e.g., 
the dynamical simulations of Rasio & Shapiro 1992, 1994, 
[ 1999; Shibata & Uryu 2000, 2002; Faber, Rasio & Manor 
■ 2001; Oechlin, Rosswog & Thielemann, 2002; Shibata, 
\ Taniguchi & Uryti 2003; Faber, Grandclement & Rasio 
. 2003; see also the review of Baumgarte and Shapiro 2003). 
It is likely that differential rotation will play an important 
role in the dynamical stability of these remnants, since 
> it can be very effective in increasing the their maximum 
allowed mass (Baumgarte, Shapiro & Shibata 2000, here- 
after BSS; Lyford, Baumgarte & Shapiro 2003, hereafter 
LBS). 

Most neutron stars in binaries have individual gravi- 
tational masses close to 1.4 Mq (compare Table 2 be- 
low). Furthermore, most recent realistic nuclear equations 
of state predict a maximum allowed mass for nonrotating 
neutron stars in the range of about 1.7 - 2.3 Mq (com- 
pare Table 1 below). Taken together, these facts seem to 
suggest that the coalescence of binary neutron stars would 
lead to prompt collapse to a black hole. However, rota- 
tion, and especially differential rotation, can increase the 
maximum allowed mass significantly^. 

The maximum allowed mass of uniformly rotating stars 
is limited by the spin rate at which the fluid at the equator 
moves on a geodesic (the Kepler limit); any further speed- 
up would lead to mass shedding. Uniform rotation can 
therefore increase the maximum allowed mass by about 20 



% at most for very stiff equations of state (Cook, Shapiro 
& Teukolsky, 1992, 1994a, 1994b, hereafter CSTl, CST2 
and CST3 respectively; see also Table 1 below), which is 
not sufficient to stabilize remnants of binary neutron star 
coalescence. Uniformly rotating neutron stars with rest 
masses exceeding the maximum allowed rest mass for non- 
rotating stars (for the same equation of state) are referred 
to as supramassive neutron stars. 

Differential rotation, however, is much more effective in 
increasing the maximum allowed mass. Unlike a uniformly 
rotating star, the rotation rate at the core of a differen- 
tially rotating star is not restricted to the maximum rota- 
tion rate at the equator, so that the core can be supported 
by rapid rotation without the equator having to exceed 
the Kepler limit. This effect was demonstrated in Newto- 
nian gravitation by Ostriker, Bodenheimer & Lynden-Bell 
(1966) for white dwarfs, and in general relativity by BSS 
for n — 1 polytropes. BSS also showed by way of illus- 
tration that stars with about 60% more mass than the 
maximum allowed mass of the corresonding nonrotating 
star can be dynamically stable against both radial and 
nonaxisymmetric modes. BSS refer to differentially rotat- 
ing equilibrium configurations with rest masses exceeding 
the maximum rest mass of a uniformly rotating star as 
hypermassive neutron stars. LBS generalized the equilib- 
rium results of BSS to other polytropic indices in the range 
0.5 <n< 2.9, and found that the largest relative increases 
in the maximum allowed mass can be found for polytropic 
indices close to n = 1. It is therefore not surprising that 



^ For coalescence from the innermost stable circular orbit, shear and shock-induced thermal pressure, which can also increase the maximum 
allowed mass, is found to have a smaller effect. 
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the hypermassive binary neutron star remnants formed in 
the n = 1 simulations of Shibata, Taniguchi & UryQ (2003) 
do not collpase promptly to black holes unless they ex- 
ceed the maximum nonrotating mass by more than about 
70%. The merged remnants will be dynamically stable 
on a dynamical timescale. Viscous damping and magnetic 
braking of differential rotation will likely occur on a sec- 
ular timescale (which is much greater than the dynami- 
cal timescale) , leading to a delayed collapse and a delayed 
burst of gravitational radiation (BSS: Shapiro 2000; Cook, 
Shapiro & Stevens 2003; Liu & Shapiro 2003). 

In this paper, we generalize the results of BSS and LBS 
to realistic equations of state. We construct differentially 
rotating neutron stars for a sample of six cold nuclear equa- 
tions of state, and tabulate their maximum allowed masses 
as a function of the degree of differential rotation. In ad- 
dition, we identify an "effective" polytropic index for each 
of these equations of state, and construct differentially ro- 
tating polytropes for these indices. We consistently find 
larger mass increases for the polytropes than for the nu- 
clear equations of state, which is due to a drop in the 
maximum density for a rotating star with a large degree 
of differential rotation. However, even for nuclear equa- 
tions of state we find increases in the maximum allowed 
rest mass on the order of about 50%, so that binary neu- 
tron star coalescence should often result in hypermassive 
neutron stars. We explicitly predict the fate of the six dou- 
ble neutron star binaries that have been identified to date, 
using the results of our equilibrium model calculations. 

The paper is organized as follows. In Section 2 we briefly 
summarize the numerical method and discuss the nuclear 
equations of state adopted in our calculations. In Section 
3 we present numerical results and compare our findings 
for nuclear equations of state with those for polytropes. 
We summarize our findings and discuss consequences for 
the six observed binary neutron star systems in Section 
4. In Appendix A we tabulate our numerical results. We 
adopt gravitational units and set G = c = 1. 

2. CONSTRUCTING NUMERICAL MODELS 

2.1. Equilibrium Models of Differentially Rotating Stars 

As in BSS and LBS, we use a modified version of the 
numerical code of GST to construct equilibrium models of 
differentially rotating neutron stars in general relativity. 
The code is based on similar numerical methods developed 
by Hachisu (1986) and Komatsu, Eriguchi, & Hachisu 
(1989), and we refer to GSTl for details. 

Gonstructing differentially rotating neutron star models 
requires choosing a rotation law F{Q) — u^u^, where u* 
and fi0 are components of the four- velocity u" and Q is the 
angular velocity. We follow GSTl and assume a rotation 
law F{Cl) = A^(Clc — ^); where the parameter A has units 
of length and where 17c is the central value of the angular 
velocity. Expressing m* and in terms of fl and metric 
potentials yields equation (42) of GSTl, or 



1 i-2f2 sin^ 6* 
in the Newtonian limit. Here we have rescaled A and 
r in terms of the equatorial radius Re'. A = A/ R^ and 
f = r/Re- The parameter ^4 is a measure of the degree of 
differential rotation and determines the length scale over 



which n changes. Since uniform rotation is recovered in 
the limit A~^ — > 0, it is convenient to parametrize se- 
quences by A~^. In the Newtonian limit the ratio be- 
tween the central and equatorial angular velocities flc/^e 
is related to A by Qc/^e = 1 + A~^, but for relativistic 
configurations this relation holds only approximately. 

We adopt this particular rotation law for convenience 
and for easy comparison with many other authors who 
have assumed the same law. However, as we pointed out 
in LBS, this rotation law approximates reasonably well 
the angular velocity profile of binary neutron star rem- 
nants in the fully relativistic dynamical simulations of Shi- 
bata & Uryu (2000, 2002) and the post-Newtonian simu- 
lations of Faber, Rasio, & Manor (2001; see also Faber, 
Grandclement & Rasio 2003). This suggests that the 
above rotation law may provide a reasonable parametriza- 
tion of differential rotation profiles that one might expect 
to find in binary neutron star merger remnants. 

We modified the numerical algorithm of GSTl by fixing 
the maximum interior density instead of the central den- 
sity for each model. This change allows us to construct 
higher mass models in some cases, since the central den- 
sity does not always coincide with the maximum density 
and hence may not specify a model uniquely. 

For each equation of state and a given value of A we 
construct a sequence of models for each value of the max- 
imum density by starting with a static, spherically sym- 
metric star and then decreasing the ratio of the polar to 
equatorial radius, TZ = Rp/Re, in decrements of 0.015. 
This sequence ends when we reach mass shedding (for 
large values of A) or when the code fails to converge (in- 
dicating the termination of equilibrium solutions) or when 
TZ = (beyond which the star would become a toroid). 
For each one of these sequences, the maximum achieved 
mass is recorded. We repeat this procedure for different 
values of the maximum density, which yields the maxi- 
mum achieved mass as a function of maximum density. 
The maximum of this curve is the maximum allowed mass 
for this particular equation of state and the chosen value 
of A. In the Appendix we tabulate our numerical results 
for the equations of state described in Section 2.2, and for 
A^^ = 0,0.3,0.5,0.7 and 1.0. Glearly, our maximum al- 
lowed masses are lower limits in the sense that even higher 
mass models may exist for other values of A or different 
differential rotation laws. 

2.2. Equations of State 

In this paper we adopt six cold, nuclear equations of 
state (EOSs) which are listed, together with numerical re- 
sults, in Table 1. The first five. A, D, L, UT and FPS are 
adopted from CST3 (who in turn adopted the labeling of 
Friedman, Ipser, & Parker (1986) for the first three EOSs). 
The last EOS, APR, was adapted from Akmal, Pandhari- 
pande & Ravenhall (1998) by Ravenhall and Pandhari- 
pande (with the smooth matching between low and high 
density provided by G. Gook). In Figure 1 we plot the 
pressure P as a function of the rest-mass density po for 
these EOSs. 

EOS A (Pandharipande 1971) models the interaction of 
neutrons at high densities with a Reid soft-core potential. 
EOS D is model V of Bethe & Johnson (1974). In EOS 
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Fig. 1. — Pressure as a function of rest-mass density for the six nuclear equations of state adopted in this paper (see text and Table 1 for 
details). 



L the nucleon interaction is modeled in terms of a mean 
scalar field (Pandharipande, & Smith 1975). Both EOSs 
UT (Wiringa, Fiks, & Fabrocini 1988) and FPS (Lorenz, 
Ravenhall, & Pethick, 1993) are modern versions of an ear- 
lier EOS proposed by Friedman & Pandharipande (1981), 
which employs both two-body (U14) and three-body nu- 
cleon interactions (TNI). EOS UT improves the treatment 
of matter at high densities, while FPS describes the inter- 
actions in terms of a Skyrme model. EOS APR (Akmal, 
Pandharipande, & Ravenhall 1998) adopts a modern two- 
nucleon interaction (A18) together with boost corrections, 
as well as the UIX thrcc-miclcon potential. 

The different descriptions of nucleon interactions af- 
fect the EOS only at high densities. At low densi- 



ties (po ^ lO^g cm^"^), the EOSs employ the Feynman, 
Metropohs, & Teller (1949) EOS, joining onto the Baym, 
Pethick, & Sutherland (1971) EOS up to neutron drip at 
po ~ 4 X lO^^g cm~^. Above neutron drip, EOSs A, D 
join onto the Baym, Bethe, & Pethick (1971) EOS, while 
EOSs L and UT join onto the Negele & Vauthcrin (1973) 
EOS. EOS APR joins onto the FPS EOS below a number 
density of 0.1 fm'^ (po = 1-26 x lO^'^g cm'^), 

All EOSs are read into our numerical code in tabular 
form, listing rest-density po, the total energy density e, 
and the pressure P at discreet points. Intermediate values 
of these quantities are computed by interpolation. Further 
details on the numerical implementation of these EOSs, in- 
cluding the fitting and tabulation, can be found in CST3. 



Table 1 

Maximum Mass Configurations 



EOS 






AfUNI 




(5Mo/Afo)"^i 


^DIF 


A//o™F 


(5Mo/A^o)°^^ 




Ab 


1.66 


1.92 


1.95 


2.24 


0.17 


2.62 


2.90 


0.51 


0.7 


D " 


1.65 


1.89 


1.95 


2.22 


0.18 


2.71 


3.02 


0.60 


0.7 




2.70 


3.23 


3.27 


3.87 


0.20 


4.45 


5.04 


0.56 


0.7 


UT " 


1.84 


2.17 


2.19 


2.55 


0.18 


2.89 


3.25 


0.50 


0.7 


FPS f 


1.80 


2.10 


2.12 


2.45 


0.17 


2.69 


3.08 


0.46 


0.5 


APR« 


2.20 


2.67 


2.46 


3.10 


0.16 


2.95 


3.50 


0.31 


0.3 



"^AU masses are in units of Solar mass Mq. 

''Raid soft core; Pandharipande (1971). 

=Model V; Bethe and Johnson (1974). 

"^Mean field; Pandharipande and Smith (1975). 

"VYIA + TNI; Wiringa, Fiks, and Fabrocini (1988). 

■''UV14 -I- TNI; Lorenz, Ravenhall, and Pethick (1993). 

9A18 +5v + UIX*; Akmal, Pandharipande, and Ravenhall (1998). 
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3. RESULTS 

Our numerical results are summarized in Table 1, 
where we list, for each equation of state, the nonrotating 
(Tolman-Oppenheimer-Volkoff, TOV) maximum masses 
(both gravitational mass M'^°^ and rest mass Mj'^^), 
the maximiim masses for uniform rotation, and the max- 
imum masses that we found in our survey of differential 
rotation. For convenience we also list fractional mass in- 
creases with respect to the TOV maximum mass, as well 
as the parameter for which the maximum mass was 
encountered. More detailed results, with results for all 
values of considered in this paper, are tabulated in 
Appendix A. 

The fractional mass increases are typically in the order 
of 50%. These increases are still large enough to support 
the remnant of a coalescing binary neutron star system 
in many cases (compare Section 4). However, they are 
noticeably below the increases for polytropic EOSs. LBS 
found that the fractional mass increases depend strongly 
on the polytropic index n, i.e. the stiffness of the EOS. 
For values of n between 0.75 and 1.25, corresponding to 
moderately stiff EOSs, they found relative mass increases 
exceeding 100%. It is therefore somewhat surprising that 
the nuclear EOSs considered in this paper yield signifi- 
cantly smaller numbers. 

To understand this result we construct models of differ- 
entially rotating neutron stars for polytropic EOSs 

P = Kpl+'^\ (2) 

where K a polytropic constant and n the polytropic index. 
For each nuclear EOS we identify an "effective" polytropic 
index ries as follows. We first construct, for a nuclear EOS, 
the nonrotating (TOV) maximum-mass model, and com- 
pute for this model the central concentration (ec/£)nuci- 
Here Cc is the central value of the energy density e, and we 
define the average energy density from 

where M is the gravitational mass and R the circumfer- 
ential radius. In Newtonian gravitation, this central con- 
densation would correspond to a unique value of n. In 
general relativity, however, this is no longer the case. To 
identify the effective polytropic index UeS, we therefore 
construct TOV maximum mass configurations for a large 
sample of polytropic indices n, and compute their central 
condensation (ec/e)poiy as described above (compare Ta- 
ble 1 in LBS). We then interpolate between the polytropic 
values of the central condensation (ec/e)poiy to the desired 
nuclear value (ec/e)nuci, and thereby identify the effective 
polytropic index ndi- We fix the polytropic constant K in 
(2) by setting the mass of the polytropic maximum mass 
TOV model, which scales with /f"/^, equal to its counter- 
part for the corresponding nuclear EOS. 

We consistently find that the agreement between the 
nuclear EOS models and the corresponding polytropes is 
quite good (to within 10% or so) for those models that 
have a maximum density Cmax similar to the maximum 
mass TOV star. However, for larger degrees of differen- 
tial rotation the maximum mass models shift to lower val- 
ues of Eniax as thc star become increasingly "toroidal" and 
bloated perpendicular to the rotation axis. We show this 
effect for EOS FPS in Figure 2, and we find qualitatively 



indentical results for all other EOSs (see also Figure 1 in 
BSS). 

This behavior can be understood by computing, as an 
alternative to our construction above, a "nuclear" poly- 
tropic index from the slope of the graphs in Fig. 1. Using 
the adiabatic index 

r = ^ (4) 



ainpo 

we difine a nuclear polytropic index by 

1 

^nucl — 



1 



(5) 



Clearly this index is a function of density, and its value 
reflects different interactions dominating different density 
regimes. Across phase transitions, rinuci can even change 
discontinously. We include a plot of npps as a function of 
the density e in the inserted panel in Figure 2. 

For a polytropic EOS, on the other hand, n is strictly 
independent of the density. We can therefore expect good 
agreement between polytropic and nuclear models only 
as long as ricn is close to rinuci for most of the mat- 
ter in thc star. As can be seen in the panel in Figure 
2, rieff is very close to npps for densities in the range 
(1.5 — 3) X lO^^gcm"^. This coincides quite well with the 
regime of maximum densities in equilibrium models for 
which the agreement between the polytropic and nuclear 
models is fairly good. For lower densities, npps drops sig- 
nificantly below Tiefi, and it is therefore not surprising that 
stellar models with maximum densities in this regime no 
longer agree very well. 

For all other EOSs we find very similar results. In all 
cases, the nuclear adiabatic index rinuci is very close to 
the effective polytropic index ries (as determined from the 
central condensation) for densities close the the maximum 
density of the maximum mass TOV star. Matter in this 
density regime dominates the structure of the star for uni- 
form and moderate degrees of differential rotation, leading 
to good agreement between nuclear and polytropic mod- 
els. For larger degrees of differential rotation, the star 
becomes increasingly "toroidal" and bloated, and the max- 
imum density shifts to smaller values. As a consequence, 
the star is dominated by matter at smaller densities, for 
which Unuci no longer agrees with UeS- Accordingly, the 
results for the polytropic models are no longer in good 
agreement with the nuclear models. 

4. SUMMARY AND DISCUSSION 

We have constructed models of differentially rotating 
neutron stars in general relativity for a sample of realistic 
nuclear equations of state. We find that their maximum 
rest masses are typically larger than those of non-rotating 
neutron stars by about 50%. This increase is significantly 
less than the increases for moderately stiff polytropic equa- 
tions of state (compare LBS). This deviation can be ex- 
plained in terms of a drop in the maximum density for 
stars with large degrees of differential rotation. At these 
smaller densities, the stiffness of the nuclear equation of 
state, as determined from equation (5), is different from 
that at the maximum density of maximum mass nonrotat- 
ing or uniformly rotating stars. Consequently, predictions 
based on simple polytropic models break down for suffi- 
ciently large degrees of differential rotation. 



EFFECT OF DIFFERENTIAL ROTATION 



5 




Fig. 2. — The maximum rest mass as a function of maximum density tmax, for different values of A. Wc include results for the FPS equation 
of state (solid lines), as well as for an n = 0.648 polytrope (dashed lines). The agreement is quite good for large values of emax, but less 
so for smaller values of emax (for rotating as well as non-rotating stars). This can be understood in terms of the polytropic index npps (as 
computed from (5); solid line in the small panel) which drops significantly below rieff (as computed from central condensation; dashed line). 
The raggedness of some of the lines is a consequence of the finite parameter step size in our sequences, and is a measure of the numerical 
error in our results. 




Fig. 3. Rest mass Mo as a function of gravitational mass M for the six nuclear equations of state (labeled as in Fig. 1). The vertical 
lines mark the gravitational mass of the individual neutron stars in those four binary neutron star systems for which the individual masses 
have been well established. 
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Even a mass increase of 50% may well be sufBcient 
to support the remnant of binary neutron star coales- 
cence. In fact, our findings can be used to predict pos- 
sible fates for the six known double neutron star bina- 
ries as listed in Table 2. For four of those, B1524+12, 
B1913+16, B2127+11C and the newly discovered binary 
J0737-3039, the individual gravitational masses are known 
to high accuracy (Thorsett & Chakrabarty 1999; Burgay 
et. al. 2003, Lyne et. al. 2004). For the two remaining ones, 
J1518+4904 and B2303+4G, the orbital inclination is un- 
known, so that the combined gravitational mass is known 
to much higher accuracy than the individual masses. As- 
suming no mass loss during the coalescence and merger 
(Faber, Grandclcmcnt & Rasio 2003) the total rest mass 
of the binary is conserved (the total gravitational mass is 
not conserved, due to neutrino and gravitational radiation 
losses). The rest mass corresponding to a particular gravi- 
tational mass depends on the equation of state, as is shown 
in Figure 3 for the four binaries with well-established in- 
dividual masses. 

The total rest mass Mq°* of the binary is an indicator 
of the fate of the merger remnant. If the total rest mass is 
less than the maximum (nonrotating) TOV mass 
then the remnant can be supported against gravitational 
collapse even without rotation. This possibility is identi- 
fied by "NS" in Table 2. If M^°* is greater than Mj^"^ , 
but less than the maximum mass of uniformly rotating 
stars Mj^^^ . then the nicrgc;r could lead to a supramassive 
star, supported by uniform rotation (denoted by "SNS" 
in Table 2). If M^°* is greater than M^^^ but less than 
the maximum mass wc find for differentially rotating stars 
M^^^ , merger may lead to a hypermassive neutron star 
supported by differential rotation (denoted by "HNS"). 
If, finally, M($°* is greater than M^^^ , the remnant will 
promply collapse to a black hole ("BH"). Clearly, these 
limits are only approximations, and the true outcome will 
depend not only on the mass and EOS but also on the 
non-equilibrium dynamics of the merger. 

Since the rest masses of the observed binaries depend 



on the equation of state, so do the possible fates. More- 
over, since all binary neutron star systems known to date 
have similar total masses, their fate, for a given equation of 
state, is presumably similar. In Table 2, the predictions for 
four of the six binaries are identical; only the high mass bi- 
nary B1913-I-16 and the low mass binary J0737-3039 show 
any differences for some of the equations of state. 

Table 2 suggests two immediate conclusions. For all the 
recent equations of state we survey and for all binaries ex- 
cept for B1913+16 the coalescence is likely to lead to a 
neutron star. For three out of the six equations of state, 
the merger of these binaries could lead to a hypermassive 
neutron star. Unlike supramassive neutron stars, hyper- 
massive neutron stars last only a fairly short time, until 
magnetic fields or viscous dissipation drive the objects to- 
ward uniform rotation, triggering a delayed collapse to a 
black hole (see BSS). 

Secondly, Table 2 shows that the outcome of binary 
neutron star coalescence is quite sensitive to the equation 
of state. This means that the observational identification 
of a merger remnant, and its possible delayed collapse to 
a black hole, could place a tight constraint on the equa- 
tion of state. Advanced LIGO (the Laser Interferometer 
Gravitational Wave Observatory) and other new gravita- 
tional wave laser interferometers may be able to observe 
binary neutron star coalescence. Comparing the most re- 
cent equations of state, UT, FPS and APR, it seems likely 
that the coalescence will lead to a rotationally supported 
neutron star. Observation of the remnant's delayed col- 
lapse to a black hole, or its absence, would then help to 
put contraints on possible equations of state. 

It is a pleasure to thank G. Ravenhall and V. J. Pand- 
haripande for providing the data for equation of state 
APR. This paper was supported in part by NSF Grant 
PHY-0139907 at Bowdoin College and NSF Grants PHY- 
0090310 and PHY-0205155 and NASA Grant NAG 5- 
10781 at the University of Illinois at Urbana-Champaign. 



Table 2 

Mass data for double neutron star binaries 



Binary 


Grav. Masses (Mq) 


Combined Mass (Mq) 


A ^ 


D 


L 


UT 


FPS 


APR 


J1518-f4904 


1.05, 1.56 


2.62 ±0.07 


BH 


HNS 


NS 


HNS 


HNS 


SNS 


B1534-I-12 ^ 


1.339, 1.339 


2.678 ±0.012 


BH 


HNS 


NS 


HNS 


HNS 


SNS 


B1913-M6 ^ 


1.3874, 1.4411 


2.8285 ±0.0014 


BH 


BH 


NS 


HNS 


BH 


HNS 


B2127+11C " 


1.349, 1.363 


2.7122 ±0.0006 


BH 


HNS 


NS 


HNS 


HNS 


SNS 


B2303-I-46 " 


1.30, 1.34 


2.64 ±0.05 


BH 


HNS 


NS 


HNS 


HNS 


SNS 


J0737-3039 ^ 


1.250, 1.337 


2.588 ±0.003 


HNS 


HNS 


NS 


HNS 


HNS 


SNS 



"Possible outcomes of binary merger for each equation of state according to computed total rest mass: (NS) neutron star; (SNS) supramassive 
neutron star; (HNS) hypermassive neutron star; (BH) black hole. 

''For these systems the binary inclination angle has not been determined precisely, so that the combined mass is known to much higher accuracy 
than the individual neutron star masses. 

•^Thorsett & Chakrabarty (1999). Lower and upper limits for comined mass are 95% central limits. 
'*yne et. al. (2004). 



EFFECT OF DIFFERENTIAL ROTATION 



7 



APPENDIX 

NUMERICAL RESULTS FOR MAXIMUM MASSES 

We list in the Tables below values for the maximum rest masses for uniformly and differentially rotating polytropes. For 
each equation of state we tabulate the differential rotation parameter , the maximum energy density emax, the ratio 
of the central and equatorial angular velocity Qc/^e, the ratio between polar and equatorial radii Rp/Rg, the ratio of the 
(relativistic) rotational kinetic energy and the gravitational binding energy T/|II^|, the maximum rest mass M™^^, and 
the fractional rest mass increase {SMo/Mq). Wc also construct differentially rotating polytropes for effective polytropic 
indices as described in Section 3 and listed in each Table heading. Wc include in the Tables the fractional mass 
increases ((5Mo/Mo)poiy for these polytropes. 

All models are computed with 64 zones both in the radial and angular direction, and with the Legendre polynomial 
expansion truncated at ^ = 16 (see CSTl for details of the numerical implementation). The error in the ratio Rp/Rg is 
determined by our stcpsize of 0.015 in this ratio (compare Section 2.1). Wc finally note that highly toroidal configurations 
depend very sensitively on the input parameters, so that those mass increases depend on the parameters adopted in our 
restricted survey. 
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Table A3 

A: M^°^ = 1.92M0, e/€c = 0.326, ries = 0.621 





Emax/lO^'"^ g Cm-3 


/ 


Rp/Re 


T/\W\ 




(SMo/Mo) 


((5Afo/Mo)poiy 


0.0 


3.55 


1.000 


0.565 


0.120 


2.24 


0.17 


0.19 


0.3 


3.45 


1.465 


0.460 


0.180 


2.50 


0.30 


0.34 


0.5 


1.87 


2.026 


0.415 


0.227 


2.69 


0.40 


0.56 


0.7 


0.938 


2.448 


0.220 


0.293 


2.90 


0.51 


0.76 


1.0 


1.12 


3.584 


0.310 


0.235 


2.34 


0.22 


0.43 
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Table A4 

D: Mj°^ = 1.89M0, e/ec = 0.284, n^s = 0-738 





(iiiax/'lO^ ' g tur ' 






T/]V\ 






(()\U|)/'.1/o)poly 


0.0 


2.73 


1.000 


0.565 


0.110 


2.22 


0.18 


0.18 


0.3 


2.64 


1.344 


0.505 


0.153 


2.41 


0.28 


0.28 


0.5 


1.51 


1.963 


0.385 


0.235 


2.95 


0.56 


0.55 


0.7 


0.884 


2.455 


0.275 


0.279 


3.02 


0.60 


1.0 


1.0 


0.774 


3.463 


2.965 


0.224E-07 


2.97 


0.57 


0.66 



Table A5 

L: Mj°^ = 3.23M0, e/ec = 0.339, rieff = 0.586 

£^ax/10^' g Cm-3 Rp/Re T/\W\ M^'^^ / Mq (SMo/Mq) (Mfo/Mo)poiy 



0.0 


1.26 


1.000 


0.550 


0.134 


3.87 


0.20 


0.19 


0.3 


1.13 


1.505 


0.425 


0.210 


4.52 


0.40 


0.37 


0.5 


0.498 


1.970 


0.320 


0.273 


4.78 


0.48 


0.52 


0.7 


0.378 


2.502 


0.010 


0.298 


5.04 


0.56 


0.69 


1.0 


0.374 


3.464 


0.025 


0.272 


4.06 


0.26 


0.35 



Table A6 







UT: Afo^ov 


= 2.17Mq 




0.321, neff = 


0.635 






emax/lO^'"' g cm-^ 


^ j ^\ p 


Rp/Re 


Tim 


Mr^'/MQ (SMo/Mo) 


((5Mo/Mo)poiy 


0.0 
0.3 
0.5 
0.7 
1.0 


2.73 
2.66 
2.02 
0.763 
0.744 


1.000 
1.455 
2.145 
2.441 
3.468 


0.565 
0.460 
0.490 
0.220 
0.010 


0.122 
0.182 
0.195 
0.293 
0.273 


2.56 
2.86 
2.91 
3.25 
2.97 


0.18 
0.32 
0.34 
0.50 
0.37 


0.19 
0.33 
0.57 
0.81 
0.47 



Table A 7 

FPS: Mj°^ = 2.IOM0, e/ec = 0.316, n^s = 0.648 
A-^ emax/lO'' g cin-3 njfle Rp/Re T/\W\ M^'^^/Mq {SMq/Mq) ((5Mo/Mo)poiy 



0.0 


2.92 


1.000 


0.565 


0.117 


2.45 


0.17 


0.19 


0.3 


2.83 


1.440 


0.475 


0.172 


2.72 


0.29 


0.32 


0.5 


1.36 


1.942 


0.360 


0.248 


3.08 


0.46 


0.59 


0.7 


0.836 


2.446 


0.275 


0.278 


3.06 


0.45 


0.82 


1.0 


0.748 


3.469 


0.010 


0.275 


2.99 


0.42 


0.47 



Table A8 







APR: Mo™^ 


= 2.67M0 


, e/ec = 


0.377, neff = 


0.495 




i-l 


Emax/lO^^ g cm" 


/ 


Rp/Re 


T/\W\ 


ms^'^'/Mq 


((5Mo/Mo) 


((5Afo/Mo)poiy 


0.0 


2.44 


1.000 


0.580 


0.132 


3.10 


0.16 


0.20 


0.3 


1.77 


1.553 


0.445 


0.210 


3.49 


0.31 


0.42 


0.5 


1.06 


1.970 


0.370 


0.245 


3.34 


0.25 


0.76 


0.7 


0.770 


2.449 


0.280 


0.275 


3.19 


0.19 


0.46 


1.0 


0.691 


3.451 


0.010 


0.276 


3.12 


0.17 


0.81 



